Comparison of nonequilibrium processes in p+Ni and p+Au collisions at GeV energies 
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The energy and angular dependence ol double differential cross sections d^a/dQdE were measured 
for p,d,t,^''^'^}ie, ^■'''^Li, ^'^'i°Be, ^^'^B, and C produced in collisions of 1.2, 1.9, and 2.5 GeV 
protons with a Ni target. The shape of the spectra and angular distributions does almost not 
change whereas the absolute value of the cross sections increases by a factor ~ 1.7 for all ejectiles 
in this beam energy range. It was found that energy and angular dependencies of the cross sections 
cannot be reproduced by the microscopic model of intranuclear cascade with coalescence of nucleons 
and the statistical model for evaporation of particles from excited, equilibrated residual nuclei. 
The inclusion of nonequilibrium processes, described by a phenomenological model of the emission 
from fast and hot moving sources, resulting from break-up of the target nucleus by impinging 
proton, leads to very good reproduction of data. Cross sections of these processes are quite large, 
exhausting approximately half of the total production cross sections. Due to good reproduction of 
energy and angular dependencies of d^a/dQdE by model calculation it was possible to determine 
total production cross sections for all studied ejectiles. Results obtained in this work point to the 
analogous reaction mechanism for proton induced reactions on Ni target as that observed previously 
for Au target in the same beam energy range. 
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I. INTRODUCTION 

The recent successful analysis [TJ |2] of the inclusive 
spectra and angular distributions of double differen- 
tial cross sections cPa/d£ldE for light charged particles 
(LCP's), i.e. particles with Z < 2, and intermediate mass 
fragments (IMF's), i.e. ejectiles heavier than ''He, pro- 
duced in proton - Au collisions at proton beam energies 
1.2, 1.9, and 2.5 GeV indicated, that the competition of 
two different processes is essential for understanding of 
the mechanism of proton - gold nucleus collisions. The 
traditional picture of the spallation reactions, namely the 
intranuclear cascade of nucleon - nucleon interactions fol- 
lowed by evaporation of particles from an equilibrated 
residuum of the cascade is not able to reproduce satisfac- 
torily the experimental cross sections. This fact leads to 
the conclusion that a significant contribution of nonequi- 
librium processes to the reaction mechanism is present. 

It was found that the difference between the data and 
theoretical cross sections from the above two-step model 
varies smoothly with the scattering angle and with the 
energy of the ejectile. It turned out, that this variation 
may be well reproduced assuming the isotropic emission 
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of nucleons and composite particles in the rest frame 
of two or three sources moving forward, i.e., along the 
beam direction. Such an effect has been interpreted in 
Refs. [1, 2_ as indication of break-up of the target nu- 
cleus by fast proton from the beam. In that physical 
picture the proton drills a cylindrical hole in the nucleus 
knocking out the small group of nucleons placed on the 
straight way of the proton through the nucleus. This 
group moves quickly in forward direction behaving as 
a fireball, which on its part emits nucleons and LCPs 
whereas the "wounded" nucleus may decay into two ex- 
cited prefragments which also serve as sources of LCPs 
and IMFs. 

The present study was performed with the aim to in- 
vestigate, whether the reaction mechanism observed and 
described in detail in Ref. [5] for Au target is also re- 
alized in collisions of protons with other target nuclei. 
The Ni target was used for this purpose because it has 
quite different properties than the Au target. The Ni nu- 
cleus is more than three times lighter, its N/Z ratio (~ 
1.1) is approximately 1.4 times smaller from that of the 
Au nucleus (N/Z ^ 1.5), and its binding energy (^ 8.8 
MeV/nucleon) is significantly larger than that of the Au 
nucleus (~ 7.9 MeV/nucleon). The appearance of the 
same reaction mechanism for such different target nuclei 
might suggest that this phenomenon is of a general char- 
acter, i.e., it occurs for all target nuclei in the studied 
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energy range. 

To facilitate the comparison of the results from the 
present study of the reactions in p+Ni system with re- 
sults of previous investigations concerning p+Au system 
[5], the present paper is organized in similar way as ref- 
erence [Ij. Experimental data are discussed in the next 
section, the theoretical analysis is described in the third 
section starting from IMF data and followed by analysis 
of LCPs cross sections, discussion of results are presented 
in the fourth and summary with conclusions in the fifth 
section. 



II. EXPERIMENTAL DATA 

The experiment was performed with the selfsupport- 
ing Ni target of the thickness about of 150 /ig/cm^, irra- 
diated by internal proton beam of COSY (COoler SYn- 
chrotron) of the Jiilich Research Center. The experimen- 
tal setup and procedure of data taking were in details 
described in Refs. [P and The beam was operated in 
so called supercycle mode to assure identical experimen- 
tal conditions for all three studied proton energies - 1.2, 
1.9 and 2.5 GeV, i.e., the same setup, electronics, the 
target thickness and its position. In this mode several 
cycles were alternated for each requested beam energy, 
consisting of protons injection from JULIC cyclotron to 
COSY ring, their acceleration with the beam circulating 
in the ring below the target, and irradiating the target 
by slow movement of the beam in the upward direction. 
The using of supercycle mode minimizes systematic ef- 
fects which might distort the studied energy dependence 
of the cross sections. 

Double differential cross sections (Pa/d^ldE were mea- 
sured at seven scattering angles; 16°, 20°, 35°, 50°, 65°, 
80°, and 100° as a function of energy of ejectiles for the 
following isotopes ^'^'^R, ^'^'^Re, '^''^'^Li, 7'°'i°Be, and 
loaig^ The carbon ejectiles were only charge identified. 

The absolute normalization of the cross sections was 
achieved by comparing the total production cross sections 
of ^Be particles, obtained by angle and energy integration 
of measured differential cross sections, with values of ^Be 
total production cross sections published in the compila- 
tion of Bubak et al. jj. It turned out that present values 
of total production cross sections at proton beam energy 
1.2 GeV are in perfect agreement with those obtained by 
Herbach et al. [5]. The same quality agreement was ob- 
served for total production cross sections of the present 
work at 2.5 GeV and Li, Be, and B cross sections mea- 
sured by Raisbeck et al. [5] on Ni target with 3.0 GeV 
proton beam. 

Typical spectra of isotopically identified ejectiles ob- 
tained in the present experiment are shown in Fig. [l] 
All spectra are smooth and do not change their shapes 
with increasing beam energy, however, the magnitude of 
the cross sections increases slightly with the energy. 
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FIG. 1: Typical spectra of "'He, ^Li, ^Be, and ^^B ejectiles 
(upper left, upper right, lower left, and lower right parts of the 
figure, respectively) measured at 35° for three energies of the 
proton beam; 1.2, 1.9, and 2.5 GeV, impinging on to the Ni 
target. Open circles represent the lowest energy, full squares 
- the intermediate energy, and the open triangles show the 
data for the highest energy. 



III. 



THEORETICAL ANALYSIS 



The analysis of the present experimental data was 
performed according to the same procedure as that ap- 
plied previously to the data from proton induced reac- 
tions on the Au target in the work of Budzanowski et 
al. [2\. First, the cross sections were evaluated from 
the intranuclear cascade with inclusion of possibility to 
coalesce the outgoing nucleons into LCPs, and with pos- 
sibility to evaporate the particles from the excited, equi- 
librated residua of intranuclear cascade. Such two-step 
model of the reaction mechanism is most frequently used 
in the literature for the description of spallation reac- 
tions at high proton energies. The INCL4.3 computer 
program of Boudard et al. [7] has been used for calcu- 
lations of intranuclear cascade and GEM2 program of S. 
Furihata [8, 9^ has been applied to evaluate evaporation 
cross sections. Since the data were generally underes- 
timated by two-step model a phenomenological analy- 
sis was performed in the following. This analysis was 
based on the assumption that additional processes exist 
besides the mechanisms described above. They were pa- 
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FIG. 2: Typical spectra of selected lithium, beryllium, and 
boron isotopes from p+Ni collisions measured at 35°, 80°, 
and 100° (left, middle, and right columns, respectively) for 
2.5 GeV proton beam impinging on to the Ni target. The 
detected particles are listed in the central panel of each row 
of pictures. Open circles represent the experimental data, 
and solid lines correspond to intranuclear cascade followed by 
evaporation of particles, respectively. 



FIG. 3: Typical spectra of lithium, beryllium, and boron ejec- 
tiles from p+Ni collisions measured at 35°, 50°, and 100° (left, 
middle, and right columns, respectively) for 2.5 GeV proton 
beam impinging on to the Ni target. The detected particles 
are listed in the central panel of each row of pictures. Open 
circles represent the experimental data, dashed, dot - dashed, 
and solid lines correspond to slow emitting source, fast emit- 
ting source and the sum of both contributions, respectively. 



rameterized by incoherent sum of isotropic emission of 
particles from highly excited sources moving in forward 
direction, i.e., along to the beam. Each of the sources 
has a Maxwellian distribution of the energy E available 
for the two body decay resulting in emission of the de- 
tected particles; cPa/dEdn - Vi^exp(-£;/T). The ve- 
locity of the source - (3 (in units of speed of light), its 
temperature - T (in MeV), and the contribution to the 
total production cross section - a are treated as free pa- 
rameters. The presence of the Coulomb barrier, which 
hinders emission of low energy particles, was taken into 
account multiplying the Maxwellian energy distribution 
by a smooth function P(E) corresponding to the trans- 
mission probability through the barrier. Two param- 
eters were introduced to characterize the properties of 
the Coulomb barrier: fc-parameter, i.e., height of the 



Coulomb barrier in units of the height of the barrier 
B of two charged, touching spheres of radius 1.44 A^^^; 

B = Zi Z2 eVl.44 (A}/^ + aI^^) MeV, and the ratio B/d 
, where d is a diffuseness of the transmission function 
through the barrier: P{E) = {I + exp{{E - kB)/d)))''^. 
Details of this procedure, as well as interpretation of pa- 
rameters of the model can be found in Appendix of Ref. 

HI- 



A. Intermediate mass fragments 

The shape of the spectra of IMFs is almost indepen- 
dent of the beam energy. Thus, only the cross sections 
measured at 2.5 GeV beam energy have been used for 
illustration of the quality of data description. On the 
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contrary, the shape of the spectra changes in regular way 
with variation of the detection angle, namely the spectra 
become more steep when the scattering angle increases. 
This is especially well visible when the data are compared 
at angles differing strongly as, e.g., 35° and 100°. Exper- 
imental spectra measured at three selected angles - 35°, 
80°, and 100° - are shown on Fig. [2] for most abundant 
isotopes of lithium, beryllium and boron as open circles 
together with the evaporation model predictions (solid 
lines) . Fluctuations of the theoretical spectra are due to 
the Monte Carlo method of model calculations, i.e. due 
to limited statistics of generated events. As can be seen, 
the theoretical spectra are more steep than the experi- 
mental ones and the absolute values of theoretical cross 
sections are several times smaller than the data. The for- 
mer effect is most pronounced for ^'^Li and ^Be cross sec- 
tions whereas the latter is present for all detected IMFs. 
Furthermore, the evaporation model does not reproduce 
the tendency of shape variation with the angle. 

All these facts indicate that an important contribution 
of another reaction mechanism must be added to evapo- 
ration cross sections to assure good reproduction of the 
data. This contribution is comparable in magnitude with 
the evaporation cross section for Li and becomes even 
more dominant for Be and B isotopes. Thus, the IMF's 
data have been analyzed in the frame of a phenomeno- 
logical model of two moving sources as it was done for 
p-|-Au reactions [T]. In this way the IMF's production on 
Ni target could be compared with that for the Au target. 

The parameters of two moving sources were searched 
for by fitting the two moving sources cross sections to ex- 
perimental data, which consisted of energy spectra mea- 
sured at seven angles: 16°, 20°, 35°, 50°, 65°, 80°, and 
100°. To decrease the number of parameters it was as- 
sumed that the velocity of the slow source emitting IMF's 
is equal to the velocity of the heavy residuum from in- 
tranuclear cascade, i.e., (3i— 0.005. The mean values of 
this velocity was found in calculations of intranuclear cas- 
cade to be equal to 0.0051c, 0.0049c, and 0.0047c for 
1.2, 1.9, and 2.5 GeV beam energy, respectively. It was 
checked that the modification of this parameter by 30% 
causes changes of other parameters smaller than their 
errors estimated by fitting computer program. In eval- 
uation of fc-parameter it was assumed that B - value is 
defined as the Coulomb barrier between the emitted par- 
ticles and the target nucleus. The results of the fit are 
not very sensitive to the value of the fc-parameter because 
the experimental low energy limit of the spectra is above 
the position of the Coulomb barrier for most of IMFs . 
Thus fixed values of fci — 0.75 and fc2 — 0.3 were used. 
The B/d ratio was arbitrarily assumed to be equal to 5.5. 

Values of the best fit parameters are listed in Table 
|T] The errors of the parameters are also given when the 
program searching for the best fit could estimate them. 
However, sometimes program was not able to estimate 
the errors, especially when strong ambiguities of param- 
eters were present. In such a case the values of the param- 
eters are quoted without estimation of errors. Sometimes 



it was useful to fix the values of parameters to avoid nu- 
merical problems leading to ambiguities. Such values are 
quoted in the table in parentheses. 

A very good description of the spectra of all IMF's has 
been obtained as can be judged from the inspection of 
Fig. [3] and values of the given in the Table |l] which 
usually vary between 1 and 2. 

It is obvious from comparison of Figs. [2]and|3]that the 
spectra of particles emitted from the slow source are very 
similar to the evaporation spectra of particles from the 
residuum of intranuclear cascade, however, it is not the 
case for the spectra of the fast source. They have much 
smaller slope and are dependent on the scattering angle. 
It is evident that their presence is necessary to reproduce 
the high energy part of experimental spectra. 

The cross sections of both emitting sources increase 
with beam energy. This is illustrated on Fig|4] where 
ratios of the production cross sections for beam energy 
1.2 and 1.9 GeV to the cross sections for 2.5 GeV are 
presented. It turned out that the increase of the cross 
sections with the beam energy is the same, in the lim- 
its of errors, for contributions emerging from both emit- 
ting sources. The cross sections increase approximately 
by factor 1.7 when beam energy increases from 1.2 GeV 
to 2.5 GeV (cf. numbers depicted in Fig. 4]). Since 
the energy dependence is very similar for both emitting 
sources, the relative contributions of the sources remain 
constant in the studied energy range and, moreover, both 
are equal in the limits of errors: The averaged over IMFs 
ratio of cr2/(o'i + <^2) is equal to 0.50(14), 0.46(11), and 
0.48(16) for beam energy 1.2, 1.9, and 2.5 GeV, respec- 
tively. Thus, this ratio averaged over three energies is 
equal to 0.48(7). 

The parameters of the sources, which influence the 
shape of the spectra (velocity P and temperature T) 
should not change with the beam energy because the 
shape of experimental spectra is independent of the beam 
energy (cf. Fig. [ij. However, there is a distinct differ- 
ence for each energy between values of the parameters 
characterizing the slow source and the fast source. The 
velocity of the fast source f32 , averaged over IMFs is equal 
to 0.027(3), 0.027(2), and 0.023(2) for beam energy 1.2, 
1.9, and 2.5 GeV, respectively. These values are about 
five times larger than fixed velocity of the slow source 
(/3i — 0.005). The temperature parameter T2 of the fast 
source is about two times larger than the temperature 
parameter of the slow source Ti. Its values (averaged 
over IMFs) are equal to 16.8(1.2), 17.4(6), and 17.8(8) 
MeV for the fast source and 7.6(3), 8.4(6), and 8.1(5) 
MeV for the slow source at three studied beam energies 
1.2, 1.9, and 2.5 GeV, respectively. 



B. Light charged particles 

All experimental spectra of LCPs from p+Ni collisions 
contain large high energy component (cf. Figs [T| [5] and 
|6|, which cannot be reproduced by evaporation of parti- 
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TABLE L Parameters of two moving sources fitted to p+Ni data for isotopically identified IMF's and for *He: Pi, Ti, and at 
correspond to source velocity, its apparent temperature, and total production cross section, respectively. The sum a = ai + a2 
is also listed. The left part of the Table (parameters with indices "1") corresponds to the slow moving source, and the right 
part contains values of parameters for the fast moving source. The upper row for each ejectile corresponds to beam energy 1.2 
GeV, the row in the middle to 1.9 GeV, and the lowest one to the energy 2.5 GeV. 



Slow source 


Fast source 


X 

mb 


Ejectile Ti/MeV ai/mb 


ft r2/MeV (72 /mb 


■*He 7.0(2) 244(6) 
7.3(2) 269(6) 
7.9(2) 283(8) 


0.040(3) 18.1(5) 76(6) 
0.036(2) 19.4(5) 94(6) 
0.035(2) 20.1(5) 101(7) 


320(9) 7.3 
363(9) 4.7 
384(11) 4.4 


■^Li 9.1(4) 8.3(4) 
10.4(4) 11.5(6) 
9.4(5) 11.6(8) 


0.035(3) 18.6(5) 4.1(5) 
0.037(3) 19.8(5) 4.5(6) 
0.026(3) 20.5(6) 8.0(9) 


12.4(7) 1.5 
16.0(9) 1.4 
19.6(1.2) 1.3 


''Li 8.1(7) 4.9(7) 
9.6(7) 8.0(1.0) 
9.6(1.0) 5.2(1.3) 


0.022(2) 14.7(4) 6.6(9) 
0.025(3) 15.9(5) 7.0(1.2) 
0.018(2) 16.0(5) 12.3(1.8) 


11.5(1.2) 1.3 
15.0(1.6) 1.3 
17.5(2.2) 1.4 


'^Li [8.0] [0.2] 
9.4(3.8) 0.6(5) 
8.0(1.8) 1.3(4) 


0.040(8) 14.4(2.0) 1.2(5) 
0.032(6) 17.2(1.1) 3.7(1.0) 
0.029(5) 18.0(1.0) 6.4(1.5) 


1.4(5) 1.4 
4.3(1.1) 1.0 
7.7(1.6) 1.0 


'Be 8.7(1.3) 2.8(5) 
9.5(9) 4.9(6) 
11.3(7) 6.9(6) 


0.025(3) 16.8(7) 3.6(7) 
0.025(3) 19.2(9) 3.9(8) 
0.032(5) 21.4(1.3) 2.6(7) 


6.4(9) 1.4 
8.8(1.0) 1.0 
9.5(9) 1.0 


''Be [8.6] 1.1(2) 
8.3(2.0) 2.4(7) 


[0.023] 12.1(8) 1.3(2) 
0.019(7) 16.9(2.4) 1.8(1.1) 


2.4(3) 1.1 

o fif c\\ no 
O.Dl^yJ U.o 

4.2(1.3) 1.0 


'"Be 6.1(1.9) 1.0(4) 
7.9(1.7) 1.0(4) 
6.1(1.9) 1.7(6) 


0.030(16) 27.8(9.0) 0.4(2) 
[0.023] 17.6(3.1) 0.7(2) 
0.024(11) 23.0(9) 0.8(4) 


1.4(5) 1.8 
1.7(5) 1.8 
2.5(7) 1.3 


'"B [6.0] 1.6(1.2) 
[6.0] 3.3(1.2) 
[6.0] 3.2(1.4) 


0.018(4) 15.5(3.1) 1.7(6) 
[0.023] 16.4(1.5) 2.4(4) 
[0.023] 17.7(1.8) 2.7(4) 


3.3(1.4) 1.7 
5.7(1.3) 1.8 
5.9(1.5) 1.9 


''B [6.0] 3.1(9) 
[6.0] 5.2(1.2) 
[6.0] 5.9(3.6) 


[0.023] 13.4(1.2) 1.8(3) 
[0.023] 17.8(2.2) 2.3(3) 
0.016(4) 14.8(3.2) 4.1(1.8) 


4.9(1.0) 1.3 
7.5(1.3) 1.5 
10(4) 1.7 



cles from the equilibrated remnant of the intranuclear 
cascade. Thus, the nonequilibrium emission of LCPs 
must play an important role. The coalescence of nucleons 
escaping from the target nucleus together with nucleons 
taking part in the intranuclear cascade may lead to such 
an emission. Indeed, Boudard et al. [7 and Letourneau 
et al. [TD] have shown that the microscopic calculation of 
coalescence occuring when the intranuclear cascade pro- 
ceeds is able to reproduce a large part of the observed 
effect. However, it was found that the improvement of 
the description of LCPs spectra deteriorates simultane- 
ously the proton spectra because increasing of the pro- 
duction of composite particles occurs on the account of 
decreasing emission of the nucleons. This contradiction 
led the present authors to search for another nonequilib- 
rium process, which could be responsible for the observed 
enhancement of the high energy part of the LCP spec- 
tra. It was proposed P] , that the emission from a fireball, 
i.e., a fast and hot group of nucleons consisting of tar- 
get nucleons lying on the straight way of the bombarding 
proton through the target nucleus [11], can account for 
the missing nonequilibrium component of the LCP cross 
sections. A sum of the coalescence of nucleons and the 



emission of LCPs from the fireball was found to be cru- 
cial for the very good description of the data for proton 
induced reactions on Au target at three beam energies: 
1.2, 1.9, and 2.5 GeV [5]. Furthermore, the emission of 
the fireball, which introduces a strong rearrangement of 
the target nucleus can lead to a break-up of the rem- 
nant of the target and thus to an appearing of two mov- 
ing sources also emitting intermediate mass fragments 
and LCps. Therefore, this hypothesis explains simulta- 
neously the presence of the nonequilibrium emission for 
IMFs which was discussed above. 

Since the fireball contains only several nucleons its con- 
tribution is present only for LCPs. On the contrary, the 
fast and slow excited prefragments of the target may emit 
IMFs as well as LCPs. In the present analysis their con- 
tribution to spectra of LCPs has been, however, neglected 
because it was estimated (by extrapolation of parameters 
found for IMFs to lighter ejectiles) to be much smaller 
than contributions of other reaction mechanisms. Mag- 
nitude of cross sections for emission of LCPs from two 
sources - products of break-up - were estimated to be 
~ 10% of the total cross sections, i.e., to be of order of 
errors of the fitting procedure. 
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FIG. 5: Typical spectra of protons, deuterons, and tri- 
tons (upper, middle, and lower rows of the figure, respec- 
tively) measured at 20°, 65°, and 100° (left, middle, and right 
columns of the figure, respectively) for 2.5 GeV proton beam 
impinging on to the Ni target. Open circles represent the 
experimental data, dashed, dot - dashed, and solid lines cor- 
respond to the two-step model, the emission from the fireball 
and the sum of both contributions, respectively. 



FIG. 4: Ratio of production cross section determined at 1.2 
(open circles) and 1.9 GeV (full dots) to the cross sections 
found at 2.5 GeV beam energy. The lower panel presents 
cross sections for slow sources, panel in the middle for fast 
sources, and upper panel for sum of both contributions. The 
horizontal lines depict values of the ratios averaged over IMFs; 
the solid lines for the 1.9 GeV , and the dashed lines for 1.2 
GeV. 



In the present study the same procedure of the descrip- 
tion of LCP spectra as that in ref. has been applied. 
The INCL4.3 computer program |7] has been used for 
the description of the intranuclear cascade of nucleon- 
nucleon collisions with inclusion of coalescence of nucle- 
ons, whereas the GEM2 computer program [51 [S] served 
for evaluation of evaporation of particles from heavy tar- 
get residuum remaining after the intranuclear cascade. 
The default parameter values, proposed by the authors 
of both programs, have been used, respectively. 

Since there is no explicit room for the presence of the 



fireball in the microscopic calculations performed accord- 
ing to the intranuclear cascade model, the inclusion of 
fireball emission should be accompanied by decreasing 
the contribution from direct processes simulated by in- 
tranuclear cascade and coalescence of escaping nucleons. 
Thus, the spectra of protons evaluated from intranuclear 
cascade with inclusion of coalescence and with contri- 
bution of evaporation of particles were multiplied by a 
factor F, common for all scattering angles, treated as a 
free parameter and then added to the contribution from 
the fireball emission calculated according to the formula 
of single moving source emitting isotropically the LCP's 
[T2| . The same, fixed value of the scaling factor has been 
used for further analysis of data for other LCPs. The 
parameters of the single moving source - the fireball, i.e. 
its temperature parameter - T3, velocity of the source 
- /33, total production cross section associated with this 
mechanism - (T3 were treated also as free parameters. 

Parameters (the height of the Coulomb barriers in 
units of i? - Coulomb barrier between the ejectile and the 
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TABLE II: Parameters of the fireball fitted to p+Ni data; /Js, Ts, and 0-3 correspond to fireball velocity in units of speed of 
light, its apparent temperature, and total production cross section, respectively. The upper row for each ejectile corresponds to 
beam energy 1.2 GeV, the row in the middle to 1.9 GeV, and the lowest one to the energy 2.5 GeV. Parameter F is the scaling 
factor of coalescence and evaporation contribution extracted from fit to the proton spectra. The numbers in parentheses show 
fixed values of the parameters. Note, that for the a particles contribution of two additional moving sources should be added 
with parameters given in Table |l] The columns described as F*aiNCL and F*aGEM contain total production cross sections 
due to intranuclear cascade with the coalescence and due to evaporation from the target residuum, respectively. The total 
production cross section obtained by summing of all contributions is depicted in the column denoted by a. In the case of alpha 
particles it contains also the contribution of the emission from slow and fast sources listed in Table |l] 



H/JfcJCllic 


/33 




0"3 


r 


F*aiMCL F*0-GEM 


(T 


X 






MeV 


mb 




mb 


mb 


mb 




P 


0.149(12) 38.9(2.1) 


1071(61) 


0.70(3) 


1094 


994 


3159(61) 


179 




0.156(10) 41.7(1.9) 


1222(53) 


0.70(2) 


1139 


1005 


3366(53) 


95.1 




0.163(8) 


43.2(1.5) 


1343(44) 


0.79(2) 


1286 


1123 


3752(44) 


48.8 


d 


0.105(4) 


32.5(8) 


181(5) 


[0.70] 


202 


174 


557(5) 


9.5 




0.099(3) 


33.7(6) 


234(5) 


[0.70] 


201 


196 


631(5) 


3.2 




0.100(4) 


35.8(8) 


272(7) 


[0.79] 


220 


225 


717(7) 


7.4 


t 


0.062(3) 


21.9(6) 


41.9(1.6) 


[0.70] 


40.9 


28.9 


111.7(1.6) 


1.3 




0.055(3) 


23.8(6) 


59.8(2.1) 


[0.70] 


41.3 


34.1 


135.2(2.1) 


1.4 




0.054(3) 


25.0(6) 


71.5(2.2) 


[0.79] 


45.3 


39.4 


156.2(2.2) 


1.3 


^He 


0.046(2) 


22.9(5) 


43.1(1) 


[0.70] 


31.2 


32.6 


106.6(1) 


3.4 




0.039(2) 


23.5(4) 


60.5(1.2) 


[0.70] 


31.6 


37.5 


129.6(1.2) 


3.3 




0.040(2) 


25.0(5) 


69.4(1.4) 


[0.79] 


34.5 


43.1 


147.0(1.4) 


2.9 










[0.70] 


16.8 


277 


614(9) 


7.3 










[0.70] 


16.6 


281 


661(9) 


4.7 










[0.79] 


18.1 


312 


714(11) 


4.4 



target nucleus) and parameter B/d describing difFuseness 
of the transmission function through the Coulomb bar- 
rier are fixed at arbitrarily assumed values 0.07 and 4.8, 
respectively. Values of the fitted parameters are collected 
in the Table |TT1 

The fit was performed for 7 scattering angles (16°, 20°, 
35°, 50°, 65°, 80°, and 100°). Results of the fit are pre- 
sented in Fig. 15] for protons, deuterons and tritons, and 
in Fig. |6] for ^e and '^He. Since the spectra at vari- 
ous beam energies almost do not differ in the shape, the 
comparison of theoretical cross sections with the data is 
shown only for one beam energy, namely for 2.5 GeV. 
The left column on both figures represents cross sections 
for 20^^, the column in the middle for 65", and the right 
column the data measured at 100°. The proton spectra 
are shown in the upper row of the Fig. [5] the deuteron 
and triton data in the middle and lower row, respectively. 
The ■^He data are depicted in the upper row of Fig. [6] 
whereas the ''He cross sections occupy the lower row of 
this Figure. 



It is obvious from inspection of Figs. |5] [6| that a very 
good description of the experimental data for all LCPs 
was achieved. It should be emphasized, that the values 
of the best fit parameters vary smoothly from ejectile to 
ejectile as well as from the one beam energy to another, 
thus the same mechanism seems to be responsible for the 
nonequilibrium processes for all these particles. It was, 
however, found that values of the parameters, necessary 
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FIG. 6: Same as on Fig. [5] but for ^He and ''He. The dotted 
line for *He denotes the contribution of an additional, slowly 
moving source. 



to describe the alpha particle data differ from those for 
lighter LCPs. It turned out that (i) it is necessary to use 
two emitting sources instead one fireball for good data 
reproduction, and (ii) the parameters of these sources 
have quite similar values as those for IMFs (cf. Table 
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|!]). For this reason, it may be concluded that the alpha 
particles behave rather like IMFs than as LCPs. 

Values of fireball velocity f}^ and its temperature pa- 
rameter do not change systematically with the beam 
energy and their fluctuations are so small, that it is pos- 
sible to assume that they do not change with the energy. 
Similar situation was observed for IMFs. Thus the en- 
ergy averaged values of velocities and temperature pa- 
rameters of all sources are collected in one figure - Fig. [7] 
- to allow for discussion of their dependence on the mass 
of ejectiles. 
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FIG. 7: In the lower panel of the figure the apparent temper- 
ature of the moving sources, averaged over beam energies is 
drawn as a function of the ejectile mass. Open circles and full 
dots represent values of temperature parameters T2 and Ti for 
fast and slow source, respectively. Full squares indicate tem- 
perature Ts of the fireball. The solid and dashed lines were 
fitted to the points representing the IMF's and ''He. Dash 
dotted line was fitted to points representing the LCP's. In 
the upper panel of the figure the dependence of the beam en- 
ergy averaged velocity of the sources is drawn versus mass of 
ejectiles. The symbols and lines have the same meaning as 
for the lower part of the figure with one exception: The full 
dots are not shown because the velocity of slower source was 
fixed during analysis (at velocity /3i =0.005 of heavy residuum 
of target nucleus after intranuclear cascade) and it is repre- 
sented by solid line in the figure. 



The velocities and temperature parameters are 



grouped in three distinctly separated sets corresponding 
to the slow source and Ti), to the fast source (/32 and 
T2), and to the fireball {[5^ and T3). The mass depen- 
dence - approximated by straight lines - is also different 
for each source. The slope of the mass dependence is the 
smallest for the slow source, of intermediate value for the 
fast source, and the largest for the fireball. The linear de- 
pendence of the temperature parameter on the mass of 
ejectiles is expected as a result of momentum conserva- 
tion, i.e., the recoil of the source of a given mass appears 
during the emission of ejectiles of various masses. This 
linear dependence of the temperature parameter on the 
mass of ejectile allows for the extraction of recoil cor- 
rected temperature of the source - t and the estimation 
of the mass of the source A5. If the mass of the source is 
constant - the same for all ejectiles - the determination of 
its mass and recoil corrected temperature is unambigu- 
ous. However, there are arguments that the source must 
have some distribution of masses with different average 
value for each ejectile. For example, the deuterons cannot 
be emitted by fireball consisted of only two nucleons but 
such emission may occur from the fireball built of three 
nucleons. The emission of protons, on the contrary, can 
appear both from two-nucleon and three-nucleon source, 
thus the mass of the fireball emitting protons is in average 
smaller than the mass of fireball emitting the deuterons. 
This may be a reason of strong decreasing of the fire- 
ball velocity with the mass of the ejectile as well as only 
slight decreasing of the velocity of the fast source with 
the mass of IMFs. A change of the mass of the fireball 
by one nucleon is very significant, because the fireball 
may be built of only several nucleons, whereas such a 
change for the source consisted of 20 or more nucleons 
is hardly to be observed in the mass dependence of the 
temperature parameter. 

The given above arguments show that the extraction 
of the mass of the source and its recoil corrected tem- 
perature from the ejectile mass dependence of the tem- 
perature parameter should be taken with caution and 
treated only as a crude estimation. Such an estimation is 
discussed below and extracted parameters are compared 
with those, which were found in our previous study [2] of 
reactions induced by protons on the gold target. 



IV. DISCUSSION 

The parameters of linear functions describing the de- 
pendence on the ejectile mass of the temperature param- 
eter T and velocity /3 of three sources are collected in the 
Table |III| Those parameters, obtained in the previous 
study of reactions induced by protons on the gold target 
[2], are also listed in this Table. 

It is clearly visible that all properties of the fireballs 
for both targets are identical in the limits of errors. This 
seems to be unexpected, especially as concerns masses 
of both fireballs, because of large difference between the 
mass of Ni and Au targets. However, according to the 
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TABLE III: Beam energy averaged temperature and velocity 
parameters of three sources of ejectiles for An and Ni targets. 
T denotes apparent source temperature (in MeV), r - the 
temperature parameter corrected for the recoil, As represents 
mass number of the source, and /3 its velocity in units of speed 
of light. The symbol A indicates the mass number of the 
ejectile. Parameters with index 1 correspond to slow source, 
with index 2 to fast source, and with index 3 to the fireball. 



Parameter 


Au 


Ni 




11.1(3) 


11.2(7) - 0.4(2)*A 


Tl 


11.1(3) 


11.2(7) 


Asi 


~ 165 


28(15) 


/3i 


[0.003] 


[0.005] 


T2 


30.6(4) - 1.61(45)*A 


22.5(6) - 0.8(1)*A 


T2 


30.6(4) 


22.5(6) 


As2 


19(6) 


28(4) 




0.059(5) - 0.0034(6)*A 


0.044(6) - 0.0021(7)*A 


Tz 


49.9(7) - 8.2(2.6)*A 


52.7(1.1) - 9.6(4)*A 




49.9(7) 


52.7(1.1) 


As3 


6(2) 


5.5(3) 


/33 


0.218(39) - 0.051(16)*A 


0.209(11) - 0.053(5)*A 



simple picture of the fireball, it is consisted of the nu- 
cleons which lie on his straight way through the target 
nucleus, therefore the mass of the fireball should scale 
as A^/^. It means that the ratio of fireball masses for 
Au and Ni should be equal to ~ 1.5. This ratio, ex- 
tracted from the phenomenological analysis, is equal to 
1.1(7) what means, that in the limits of errors it is in 
agreement with the assumed picture of the mechanism. 

The equality of velocities of fireballs and their temper- 
atures for both targets may be treated as consequence of 
the same momentum and energy transfer from the bom- 
barding proton to the group of nucleons forming the fire- 
ball. Such an explanation is in line with the fact of the 
same beam energies for both targets and small difference 
in the thickness of the nuclear matter placed on the way 
of bombarding proton. Therefore, this equality may be 
interpreted as the argument in favor of assumed model 
of the reaction. 

The recoil corrected temperatures of slow sources for 
Ni and Au targets are also the same. Of course, the mass 
As^ of the slow source is completely different in the case 
of Au target (~ 165) and Ni target (-- 28 ). Therefore the 
recoil correction of the temperature parameter could be 
neglected for the Au target but is visible in the ejectile 
mass dependence of the temperature parameter Ti for 
the Ni target. This difference of source masses reflects 
also on the velocity of the slow source /3i . 

The largest differences appear for the fast source. It 
should be, however, pointed out that the parameters do 
not differ in average more than ~ 50%. Taking into con- 
sideration the fact of limited accuracy of extraction of 
values of the parameters it may be claimed that the pa- 
rameters of the fast source are similar for both targets. 

The above considerations concern velocity and temper- 



TABLE IV: Averaged over IMFs ratio of cross section ai 
for energy E to the same cross section measured at proton 
beam energy 2.5 GeV on Au and Ni targets. The ai is equal 
to (Tl (slow source), (T2 (fast source) or a = ai + a2 (total 
production cross section) 



Ratio of ai to ai(2.5 GeV) 


at E/GeV 


Au 


Ni 


ai 


1.2 


0.39(3) 


0.63(8) 




1.9 


0.75(6) 


0.85(12) 


cr2 


1.2 


0.23(3) 


0.63(9) 




1.9 


0.57(6) 


0.82(11) 


ai + (72 


1.2 


0.33(2) 


0.57(6) 




1.9 


0.66(4) 


0.81(8) 



ature of three sources for Au and Ni targets and show, 
that this parameters, which are almost independent of 
the beam energy in the studied proton energy range 1.2 
- 2.5 GeV, are very similar for both targets. 

In the following the behavior of production cross sec- 
tions will be discussed. As it was shown on Fig. [4] cross 
sections for all IMFs increase in average by factor ~ 1.7 
when proton beam energy increases from 1.2 GeV to 2.5 
GeV. This is true for the total production cross sections 
as well as for the contributions of individual emitting 
sources, however, the spread (among various IMFs) of 
the ratio of given cross section to that measured at 2.5 
GeV is smaller for total production cross section than for 
individual sources. The ratios of the cross sections mea- 
sured at 1.2 and 1.9 GeV to those measured at 2.5 GeV 
are listed in Table HVl 

It is clear from examination of Table II VI that the IMF 
production cross sections measured for Au target increase 
stronger with beam energy than those for Ni target. This 
can be understood as an effect caused by difference be- 
tween threshold energies for fragmentation of both tar- 
gets. To illustrate this effect the excitation function for 
^Be production is shown on Fig. [sjas a typical example. 
It is seen that fragmentation starts at lower energies on 
Ni target than on Au target. Therefore the studied range 
of beam energy (1.2 - 2.5 GeV) corresponds for Ni target 
to the region where the production cross section starts to 
saturate, whereas for Au target this is region where the 
production cross section starts to rise quickly. Further- 
more, the leveling of the production cross section for the 
Ni target appears at lower value than that for Au target. 
Both these effects cause that the total production cross 
sections should rise more quickly for the Au target than 
for Ni target, in accordance with present observations. 

It is interesting to examine how large are the contri- 
butions of individual reaction mechanisms for emission 
of LCPs, i.e., the preequilibrium stage of the reaction 
described by intranuclear cascade and coalescence of nu- 
cleons into composite particles, the fireball emission, and 
the evaporation. As it is seen in Fig. [9] (right column) 
the contribution from intranuclear cascade and coales- 
cence is for Ni target almost equal to that from evapo- 



10 




1000 
Ep [MeV] 



10000 



FIG. 8: Energy dependence of production cross section of ^Be 
ejectiles in proton induced reactions. The lines siiow results of 
compilation of ^Be production cross sections [4], the symbols 
represent experimental data of the present experiment (open 
squares) for Ni target and the data published by Budzanowski 
et al. 2 (full squares) for An target. The solid and dashed 
lines depict the excitation functions from Ref. |4] for Au and 
Ni targets, respectively. 
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ration and exhausts about 30% of the total production 
cross section for all studied beam energies and ejectiles. 
The first of these contributions decreases several percent 
with the energy, whereas the evaporation contribution 
is almost independent of the energy. The contribution 
of the fireball is slightly larger (~ 40 %) and increases 
several percent in the studied energy range. These vari- 
ations are more pronounced for tritons and ^He than for 
protons and deuterons. 

This picture is quite different from energy behavior 
of separate reaction mechanisms of LCPs production for 
the Au target (cf. Fig. [9] left column), with exception of 
the evaporation contribution (cgbm on the figure) which 
is almost energy independent similarly as for Ni target. 
The coalescence contribution decreases with beam en- 
ergy (~ 10% for p,d,t, and - 30% for ^He) but fireball 
contribution increases (also ^ 10% for p,d,t, and even 
more for "^He). Moreover, the coalescence contribution 
for Au target is in average larger, and the fireball con- 
tribution smaller than for Ni target. Such a behavior 
of relative contributions of various processes might sug- 
gest that they depend rather on the proton beam energy 
per nucleon of the target than on the beam energy itself, 
similarly as it is in the case of the total production cross 
section (cf. Figs. |8]and|9]). 



FIG. 9: Energy dependence of various reaction mechanisms 
for protons (open squares), deuterons (full squares), tritons 
(open triangles) and ^He (full triangles). The relative con- 
tribution of the fast stage of the reaction, i.e., intranuclear 
cascade and coalescence of nucleons into LCPs ct/jvcl/o" is 
presented in the upper panel of the figure, evaporation con- 
tribution from the equilibrated residuum of the intranuclear 
cascade ctg-Bm/ct is shown in the middle panel, whereas the 
contribution of fireball emission aa/cr is depicted in the lower 
panel. 



V. SUMMARY AND CONCLUSIONS 

A new, extensive set of double differential cross sec- 
tions cP(j/d£ldE for the production of LCPs and light 
IMFs (^He, s.^.SLi^ ^.g.iogg^ iO'"B,and C) in collision of 
protons with Ni target has been measured at three beam 
energies (1.2, 1.9, and 2.5 GeV). The data were analyzed 
using two-step microscopic model - intranuclear cascade 
with possibility to take into account coalescence of nu- 
cleons into LCPs, followed by evaporation of particles 
from the equilibrated residual nuclei. Since the model 
could not satisfactorily reproduce the spectra measured 
for seven scattering angles; 16°, 20°, 35°, 50°, 65°, 80°, 
and 100°, the phenomenological analysis has been per- 
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formed assuming the emission of LCPs from the fireball, 
i.e., fast and hot source built of several nucleons moving 
in the forward direction, as well as emission of IMF from 
two slower and colder sources, which were interpreted as 
prefragments of the target appearing due to break-up of 
the target. An excellent description of all data has been 
achieved with smoothly varying values of the parameters 
from ejectile to ejectile. Due to such good reproduction of 
energy and angular dependencies of (Pa/dQdE by model 
calculation it was possible to determine total production 
cross sections for all studied ejectiles. They are listed in 
Table |T] and |IT] for IMFs and LCPs, respectively. 

It turned out that the competition of the emission from 
the fireball with the coalescence process and evaporation 
of particles from the equilibrated remnant of intranuclear 
cascade is crucial for a proper description of the data for 
LCPs. The emission of IMFs from a fast moving source 
was necessary to reproduce the high energy part of the 
experimental spectra. Moreover, the emission of IMFs 
from a slow moving source, which emits particles with 
spectra similar to those of evaporated IMFs from the 
equilibrated residua of intranuclear cascade, improved 
significantly the description of the low energy part of the 
spectra. 

In summary, a large contribution of the nonequilibrium 
processes, which are not contained in the two-step micro- 
scopic model, has been established. Properties of these 
processes are compatible with hypothesis of emission of 
ejectiles from three moving sources: The light, fast, and 
hot source - a fireball, which appears as result of knock- 
out of group of nucleons lying on the straight way of the 
impinging proton through the nucleus, and two slower 
and colder sources, which are created due to break-up of 



the target remnant after emission of the fireball. 

All the discussed above effects are very similar to those 
observed by present authors for another nuclear system, 
namely p-l- Au in the same range of beam energies (TJ [2] . 
As it was mentioned in the Introduction, the Ni target 
is significantly different from Au target. Thus the obser- 
vation of analogous reaction mechanisms in both nuclear 
systems suggests that this mechanism appears generally. 
Properties of this mechanism are in many aspects the 
same for Ni and for Au targets but some clear differ- 
ences are also present. Therefore further investigations 
should be performed for other nuclei and beam energies 
to achieve reasonable progress toward the understanding 
of the fragmentation phenomenon, and to allow for good 
description of the production cross sections of LCPs and 
IMFs. 



Acknowledgments 

The technical support of A.Heczko, W. Migdal, and N. 
Paul in preparation of experimental apparatus is greatly 
appreciated. This work was supported by the European 
Commission through European Community-Research In- 
frastructure Activity under FP6 project Hadron Physics, 
contract number RII3-CT-2004-506078). This work was 
also partially supported by the Helmholtz Association 
through funds provided to the virtual institute " Spin and 
strong QCD" (VH-VL231). One of us (M.F.) acknowl- 
edges gratefully financial support of Polish Ministry of 
Science and Higher Education (Grant No N N202 174735, 
contract number 1747/B/H03/2008/35). 



[1] A. Bubak, A. Budzanowski, D. Filgcs, F. Goldenbaum, 

A. Heczko, H. Hodde, L. Jarczyk, B. Kamys, M. Kistryn, 
St. Kistryn, St. Kliczewski, A. Kowalczyk, E. Kozik, P. 
Kulessa, H. Machner, A. Magiera, W. Migda, N. Paul, 

B. Piskor-Ignatowicz, M. Puchaa, K. Pysz, Z. Rudy, R. 
Siudak, M. Wojciechowski, and P. Wiistner, Phys. Rev. 
C 76,014618 (2007) 

[2] A. Budzanowski, M. Fidelus, D. Filges, F. Golden- 
baum, H. Hodde, L. Jarczyk, B. Kamys, M. Kistryn, St. 
Kistryn, St. Kliczewski, A. Kowalczyk, E. Kozik, P. Ku- 
lessa, H. Machner, A. Magiera, B. Piskor-Ignatowicz, K. 
Pysz, Z. Rudy, R. Siudak, and M. Wojciechowski, Phys. 
Rev. C 78, 024603 (2008) 

[3] R. Barna, V. BoUini, A. Bubak, A. Budzanowski, D. 
D.Pasquale, D. Filges, S. V. Fortsch, F. Goldenbaum, A. 
Heczko, H. Hodde, A. Raliano, L. Jarczyk, B. Kamys, J. 
Kisiel, M. Kistryn, St. Kistryn, St. Kliczewski, A. Kowal- 
czyk, P. Kulessa, H. Machner, A. Magiera, J.Majewski, 
W.Migdal, H.Ohm, N.Paul, B. Piskor-Ignatowicz, K. 
Pysz, Z. Rudy, H. Schaal, R. Siudak, E. Stephan, 
G.F.Steyn, R.Sworst, T.Thovhogi, M.Wojciechowski, 
W. Zipper, Nucl. Instr. Meth. in Phys. Research A 519, 
610 (2004) 



[4] A. Bubak, B. Kamys, M. Kistryn, B. Piskor-Ignatowicz, 
Nucl. Instr. Meth. in Phys. Research B 226, 507 (2004) 

[5] C.-M. Herbach, D. Hilscher, U. Jahnke, V.G. Tishchenko, 
J. Galin, A. Letourneau, A. Peghaire, D. Filges, F. Gold- 
enbaum, L. Piehkowski, W.U. Schroder, and J. Toke, 
Nuclear Physics A 765, 425 (2006) 

[6] G.M. Raisbeck, P. Boerstling, R. Klapisch, and T.D. 
Thomas, Physical Review C 12, 527(1975) 

[7] A. Boudard, J. Cugnon, S. Leray, and C. Volant, Nucl. 
Phys. A 740, 195 (2004) 

[8] S. Furihata, Nucl. Instr. and Meth. in Phys. Research B 
71, 251 (2000) 

[9] S. Furihata and T. Nakamura, Journal of Nuclear Science 

and Technology Supplement 2, 758 (2002) 
[10] A. Letourneau, A. Bohm, J. Galin, B. Lott, A. Peghaire, 

M. Enke, C. M. Herbach, D. Hilscher, U. Jahnke, V. 

Tishchenko, et al, Nucl. Phys. A 712, 133 (2002) 
[11] G.D. Westfall, J. Gosset, P.J. Johansen, A.M. Poskanzer, 

W.G. Meyer, H.H. Gutbrot, A. Sandoval and R.Stock, 

Phys. Rev. Lett. 37, 1202 (1976) 
[12] G. D. Westfall, R. G. Sextro, A. M. Poskanzer, A. M. 

Zebelman, G. W. Butler, and E. K. Hyde, Phys. Rev. C 

17, 1368 (1978) 



